Dark Matter in Modern Cosmology 



Sergio Colafrancesco 

ASI-ASDC, Via G.Galilei c/o ESRIN, 00040 Frascati, Italy 

Abstract. The presence of Dark Matter (DM) is required in the universe regulated by the standard 
general relativistic theory of gravitation. The nature of DM is however still elusive to any experi- 
mental search. We discuss here the process of accumulation of evidence for the presence of DM in 
the universe, the astrophysical probes for the leading DM scenarios that can be obtained through 
a multi-frequency analysis of cosmic structures on large scales, and the strategies related to the 
multi-messenger and multi-experiment astrophysical search for the nature of the DM. 
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DARK MATTER IN MODERN COSMOLOGY: GAINED 

EVIDENCE 

There is overwhelming evidence that we live in a geometrically flat (Q,q ^ 1) universe 
which is dominated - in the standard general relativistic cosmological paradigm - by 
a dark form of matter (Dark Matter - DM) and by an obscure form of energy (Dark 
Energy - DE). In such standard cosmological paradigm, DM provides a substantial 
fraction Qdm ~ 0.227 of the overall matter-energy content (the rest being provided by 
Dark Energy with Q.de ~ 0.728 with the baryonic contribution limited to « 0.045, 
see, e.g., Komatsu et al. 2010) and amounts to ~ 83% of the total mass content of the 
universe. 

How did we arrive to such a conclusion? Is there a general and definite consensus on it? 



The Dark Matter timeline 

The discovery of the presence of DM can be considered as a typical scientific rev- 
olution (Kuhn, 1970) that induced a change in the reference (cosmological) paradigm. 
However, as often occurs in a paradigm shift, there was no single discovery, but new 
concepts were developed and integrated step-by- step 1 . 

To begin, there are actually two dark matter problems, the local Dark Matter close to 
the plane of our Galaxy, and the global Dark Matter surrounding galaxies, clusters of 



In this context, it is necessary, since the beginning of this discussion, to acknowledge many reviews on 
the Dark Matter and related problems, as given by Faber & Gallagher (1979), Trimble (1987), Kormendy 
& Knapp (1987), Srednicki (1990), Turner (1991), Silk (1992), Ashman (1992), van den Bergh (2001), 
Ostriker & Steinhardt (2003), Rees (2003), Turner (2003), Colafrancesco (2006, 2007), Einasto (2004, 
2009), among others. 



galaxies and large scale structures. 

The local Dark Matter (Oort 1 932) in the Galactic disk is baryonic [faint stars or planet 
(jupiter)-like objects], since a collection of matter close to the galactic plane is possible if 
it has formed by the contraction of pre- stellar matter towards the plane accompanied by 
dissipation of the extra energy, so to conserve the flat shape of the population. It is now 
clear that the amount of local Dark Matter is low, and it depends on the mass boundary 
between luminous stars and faint invisible stars or planet-like objects (see Einasto 2009 
for a discussion and references therein). 

The global Dark Matter (to which we refer in this paper) is the dominating mass com- 
ponent in the universe; it is mainly concentrated in galaxies, clusters and superclusters 
of galaxies, and populates all other large-scale structures in the universe. 
The first evidence for such global DM was obtained by F. Zwicky (1933). From obser- 
vations of the radial velocities of eight galaxies in the Coma cluster, Zwicky found an 
unexpectedly large velocity dispersion o v = 1019 ± 360 km s _1 [we note in passing that 
Zwicky's velocity dispersion from only eight galaxies agrees well with the modern value 
o v = 1082 km s _1 , as obtained by, e.g., Colless & Dunn (1996)]. Zwicky concluded from 
these observations that, for a velocity dispersion of ~ 1000 km s _1 , the mean density of 
the Coma cluster would have to be ~ 400 times greater than that which is derived from 
luminous matter. 2 Zwicky concluded, therefore, that in order to hold galaxies together 
in the cluster, the Coma cluster must contain huge amounts of some dark (i.e. not visi- 
ble), cold matter (Zwicky used actually the words "Dunkle Kalte Materie" which might 
be regarded as the first reference to Cold Dark Matter ... even though not in the modern 
sense). 

Not until three years later it was found (Smith 1936) that the Virgo cluster also appears 
to exhibit an unexpectedly high mass. Smith made the speculation that the excess mass 
of Virgo "represents a great mass of internebular material within the cluster" (see van 
den Bergh 2001). 

Six years later than Zwicky's 1933 paper, Babcock (1939) obtained spectra of the 
Andromeda galaxy (M31) and found that in its outer regions the galaxy is rotating with 
an unexpectedly high velocity, far above the expected Keplerian velocity. He interpreted 
this result either as a high mass-to-light ratio in the periphery or as a strong dust 
absorption. Dark Matter in galaxies was also envisaged. 

One year later, Oort (1940) studied the rotation and the surface brightness of the edge-on 
SO galaxy NGC 3115. He found that "The distribution of mass in this system appears to 
bear almost no relation to that of light. " He concluded that a value M/L ~ 250 should 
be present in the outer regions of NGC 3115 (note that this value is reduced by almost 
an order of magnitude if the modern distance to this galaxy is adopted). Oort ended his 
paper by writing that "There cannot be any doubt that an extension of the measures 
of rotation to greater distances from the nucleus would be of exceptional interest. " 



Zwicky actually overestimated the mass-to-light ratio of the Coma cluster because he assumed a Hubble 
parameter Hq = 55%kms~ l Mpc~ l . His value for the overdensity of Coma should therefore be reduced 
from 400 to ~ 50 by using the actual value of the Hubble parameter. It is interesting to note that Hubble's 
prestige was so great at that time, that none of the early authors thought of reducing Hubble's constant as 
a way of lowering their mass-to-light ratios. 



However, no connection was made between the missing mass in this SO galaxy, and 
the Zwicky and Smith missing mass problem in clusters of galaxies. 

The DM problem remained as a kind of "anomaly" for roughly a quarter of century 
after the initial Zwicky's paper until when Kahn & Woltjer (1959) pointed out that M3 1 
and our Galaxy were moving towards each other, so that they must have completed most 
of a (very elongated) orbit around each other in a Hubble time. These authors found that 
if M31 and our Galaxy started to move apart ~ 15 Gyr ago, the mass of the Local Group 
had to be ^ 1.8 x 1O 12 M . Assuming that the combined mass of the Andromeda galaxy 
and the Milky Way system was 0.5 x 10 12 M Q , Kahn & Woltjer (1959) concluded that 
most of the mass of the Local Group existed in some invisible form. They opined that 
it was most likely that this missing mass was in the form of hot (with temperature of 
~ 5 x 10 5 K) gas. From a historical perspective, it is interesting to note that also Kahn 
& Woltjer did not seem to have been aware of the earlier papers by Zwicky (1933) and 
Smith (1936) on the missing mass problem in clusters of galaxies. 

In the 1960's and 1970's, Babcock's optical rotation curve of M31, and that of Rubin 
& Ford (1970), were extended to even larger radii by Roberts & Whitehurst (1975) using 
21 -cm line observations that reached a radial distance of ~ 30 kpc. These observations 
clearly showed that the rotation curve of M31 did not exhibit any Keplerian drop-off. 
From these observations Roberts & Whitehurst (1975) concluded that the mass-to-light 
ratio had to be ^ 200 in the outermost regions of M3 1 . Again, it is interesting to note 
that neither Babcock, nor Roberts & Whitehurst, cited the 1933 paper by Zwicky. In 
other words, no connection was made between the missing mass in the outer region of 
spiral galaxies and the missing mass in galaxy clusters such as Coma (Zwicky 1933) 
and Virgo (Smith 1936). Regarding the flat outer rotation curves of galaxies Roberts 
(1999, as quoted in van den Bergh 2001) recalls that this result "was, at best, received 
with skepticism in many colloquia and meeting presentations. " Nonetheless, the paper 
by Roberts & Whitehurst (1975) was important because it, together with papers on the 
stability of galactic disks (e.g. Ostriker & Peebles 1973), and on the apparent increase of 
galaxy mass with increasing radius (e.g. Ostriker, Peebles & Yahil 1973), first convinced 
the majority of astronomers that a missing mass problem existed. 

By the mid of 1970's enough research was done for the community to see that the 
missing mass "anomaly" was not going to go away, and the majority of astronomers 
had become convinced that missing mass existed in cosmologically significant amounts. 
"When... an anomaly comes to seem more than just another puzzle of normal science, 
the transition to crisis and to extraordinary science has begun" (Kuhn 1970). In other 
words, it began to be clear in the mid 1970's that a paradigm shift would be required to 
interpret the observations that seemed to support the ubiquitousness of missing (dark) 
matter in the universe. 

At the beginning of the 1980's the presence of DM was directly confirmed by many 
independent sources: the dynamics of galaxies and of stars within galaxies (Rubin et al. 
1978, 1980), the mass determinations of galaxy clusters based on gravitational lensing 
3 (see, e.g., Bartelmann & Schneider 2001 for a review and references therein), X-ray 



i It is noted, in passing, that Zwicky (1937) also pointed out that gravitational lensing might provide 
useful information on the total masses of galaxies. 



studies of clusters of galaxies (see, e.g., Forman et al. 1985), the indications of the first 
N-body simulations of large scale structure formation and of galaxies and cluster forma- 
tion (see, e.g., White et al. 1987). 

At this time it was also already clear that structures in the universe form by gravitational 
clustering, started from initially small fluctuations of the density of matter (Peebles 
1980). In order to form the presently observed structures, the amplitude of the density 
fluctuations must be at least one thousandth of the density itself at the epoch of recombi- 
nation, when the universe started to be transparent. The relic emission coming from this 
epoch was first detected in 1965 (Penzias & Wilson 1965) as an (almost) isotropic cos- 
mic microwave background (CMB). When finally the fluctuations (anisotropics) of the 
CMB were measured by the COBE satellite, first an large angular scales (Smooth et al. 
1992), and by the Boomerang experiment, also on smaller angular scales (de Bernardis 
et al. 2000), they appeared to be two orders of magnitude lower than expected from 
the density evolution of the luminous (baryonic) mass. The solution of the problem was 
suggested independently by several theorists. If we suppose that the dominating mass 
component of the universe - Dark Matter - is not made of ordinary matter but of some 
sort of non-baryonic matter, then density fluctuations can start to grow much earlier, 
and have at the time of recombination the amplitudes needed to form the structures ob- 
served today. The interaction of this non-baryonic matter with radiation is much weaker 
than that of ordinary matter, and radiation pressure does not slow the early growth of 
fluctuations. 

However, beyond the indubitable success of the DM-dominated cosmological frame- 
work, the fundamental uncertainty on the DM nature remained. 

The first suggestion for a non-baryonic DM was made referring to particles well known 
at that time to physicists, i.e. neutrinos (v) (see, e.g., Pontecorvo 1967), for which "the- 
oretical and practical considerations" push progressively some physicists to agree on a 
non zero mass neutrino (Quantum mechanics allow then the oscillation of neutrinos: a 
V e can, along its travel in the universe, become a v M and vice- versa). This can be con- 
sidered as the first input from particle physics to the search for the nature of DM. Rees 
(1977) was one of the first astrophysicists to tackle with the idea of non-baryonic nature 
for the cosmologically relevant DM, and suggested that: "There are other possibilities 
of more exotic character - for instance the idea of neutrinos with small f~ few eV) rest 
mass has been taken surprisingly seriously by some authors ". However, this scenario of 
hot DM soon led to major problems, because neutrinos move with very high velocities 
which prevents the formation of small structures as galaxies (Silk 1968). 

Thus some other hypothetical non-baryonic particles were suggested, such as axions 
(Abbott & Sikivie 1983, Preskill, Wise & Wilczek, 1983; Dine & Fischler, 1983) or 
Weakly Interacting Massive Particles (e.g. WIMPs, Goldberg 1983, Ellis et al. 1983). 
The essential character of these particles to make them good for cosmology is that they 
have much lower velocities than neutrinos. Because of this, the new version of Dark 
Matter scenario was called Cold, in contrast to neutrino-dominated Hot Dark Matter 
scenario. 

Numerical simulations of the evolution of the structure of the universe confirmed the 
formation of filamentary superclusters and voids in the Cold Dark Matter scenario of 
structure formation (e.g. Ostriker 1993 and references therein). The suggestion of the 
Cold Dark Matter has solved most problems of the new cosmological paradigm. 



One unsolved problem remained. Estimates of the total matter density (ordinary plus 
DM) yield values Q. m ~ 0.27 of the critical density. This value - not far from unity, but 
definitely smaller than unity - is neither favoured by theory nor by the data, including 
the measurements of the cosmic microwave background, the galaxy dynamics and the 
expansion rate of the universe obtained from the study of distant type la supernovae 
(SNe). To fill the matter/energy density gap between unity and the observed matter 
density it was assumed that some sort of Dark Energy (DE) exists. 4 By the early 
2000's, refined CMB anisotropy experiments (e.g., Boomerang, Maxima, DASI, and 
then WMAP) demonstrated that the CMB anisotropy spectrum is the largest detector 
for the presence of Dark Matter in their fluctuation spectrum, and by combining their 
measurements with independent cosmological distance measurement using type la SNe 
(Perlmutter 2000), a first determination of the overall matter-energy composition of 
the universe was possible (see Komatsu et al. 2010 for a recent determination of the 
cosmological parameter set). The inclusion of the DE term in the general relativistic 
cosmological scenario has filled the last gap in the modern cosmological paradigm. 



False alarms and diversionary manoeuvres 

No good detective story is complete without at least one false clue or some diversion- 
ary manoeuvres. 

Oort (1960, 1965) believed that he had found some dynamical evidence for the 
presence of missing mass in the disk of the Galaxy. If true, this would have indicated 
that some of the dark matter was dissipative in nature. However, late in his life, Jan 
Oort confessed (as reported by van den Bergh 2001) that the existence of missing mass 
in the Galactic plane was never one of his most firmly held scientific beliefs. Detailed 
observations, that have been reviewed by Tinney (1999), show that brown dwarfs cannot 
make a significant contribution to the density of the Galactic disk near the Sun. 

The presence of large amounts of matter of unknown origin has given rise to specu- 
lations on the validity of the Newton law of gravity at very large distances. One of such 
attempts is the Modified Newton Dynamics (MOND) model, suggested by Milgrom 
(1983) and Milgrom & Bekenstein (1987) (see, e.g., Sanders 1990 for a discussion). 
Indeed, MOND is able to explain spiral galaxies quite well imposing a minimum ac- 
celeration scale ciq ~ cHq, without assuming the presence of some hidden matter. If the 
MOND scale ciq is allowed to run, galaxy clusters might be explained as well. However, 
there exist several arguments which make this model unrealistic (see Einasto 2009 for 
a discussion). In addition, a full relativistic theory is needed in any case to construct a 
consistent cosmology. A tensor-vector-scalar theory (TeVeS, see Bekenstein 2004) has 
been then proposed and it is successful in reproducing MOND in the proper limit even 
for photons, and respecting also the classical tests of general relativity (GR). 
Whether a consistent cosmology can be constructed with this theory remains to be seen 
(in fact, both MOND and TEVES models are not covariant formalisms of a general 



This assumption was not new: already Einstein added to his cosmological equations a cosmological 
constant term A (corresponding to a vacuum energy). 



relativistic theory of gravitation, see Capozziello et al. 2008). 

Extended Gravity scenarios have been therefore suggested to explain the large-scale 
cosmological problem of the accelerated expansion of the universe and have been also 
worked out to describe the morphology and dynamics of very peculiar systems, like the 
bullet cluster (Brownstein & Moffat 2007). Such a scenario is today the best alternative 
to a pure particle DM scenario. 

The Astro-Particle connection 

All the available information indicates that the standard scenario for structure for- 
mation in the GR framework requires that the global Cold Dark Matter must be non- 
baryonic, its density fluctuations start to grow much earlier than perturbations in the 
baryonic matter, and have at the recombination epoch large enough amplitudes to form 
all structures seen in the universe. 

However, the actual nature of the CDM particles still remains unknown. Physicists have 
attempted to discover particles which have the properties needed to explain the struc- 
ture of the universe, but so far without success. This means that a true Astro-Particle 
connection should be developed in the search for the nature of DM. In conclusion, even 
though the direct information on the dark components of the universe (DM and DE) 
comes solely from astronomical observations, it is clear that a definite understanding of 
the nature of DM will come through the discovery and the multi-messenger study of the 
fundamental particles of which DM consists of. 



A last remark 

The discovery of DM has the general character of a typical scientific revolution 
connected with changes of paradigms, as discussed by Kuhn (1970) in his book "The 
Structure of Scientific Revolutions". There are not so many areas in modern astronomy 
and cosmology where the development of ideas can be described in these terms. 

DARK MATTER IN MODERN COSMOLOGY: THE PRESENT 

Motivations 

All the reliable indications for the presence of the dark components of the universe, 
and especially of DM, come solely from astronomical evidence. 

The main indications are: i) galaxy rotation curves; ii) dwarf galaxy mass estima- 
tors; iii) galaxy cluster mass estimators; iv) lensing reconstruction of the gravitational 
potential of galaxy clusters and large scale structures; v) the combination of global geo- 
metrical probes of the universe (e.g., CMB) and distance measurements (e.g., SNe). 
Supporting evidence comes also form Large Scale Structure simulations for the leading 
structure formation scenario. 



Dark Matter candidates 



There are five basic properties that DM candidates must have: DM must be dissipa- 
tionless, collisionless, cold, must behave like a fluid on galactic scales and above, must 
behave sufficiently classically to be confined on galactic scales (see Baltz 2004). The 
first three properties do not place any stringent constraint on the space of possibilities, 
while the last two place upper and lower (loose) bounds, respectively, on the mass of the 
particle. Such wide space of possibilities allowed theoreticians to propose many candi- 
dates for the DM particles (see Feng 2010). Among these, the most viable and widely 
considered candidates for the DM are, so far, neutralinos (the lightest particles of the 
minimal super symmetric extension of the Standard Model, MSSM, see e.g. Jungman 
et al. 1996), sterile neutrinos (the lightest right-handed neutrino, see e.g. Shaposhnikov 
2007) or even other forms of light DM (see e.g. Boyanovsky et al.2007). 
In the following, and for the sake of brevity, we will focus mainly on the astrophysi- 
cal probes related to some of these DM candidates, specifically neutralinos and sterile 
neutrinos. We will also focus our discussion, for the sake of conciseness, on the best 
astrophysical laboratories for DM indirect search: i) galaxy clusters (the largest bound 
containers of DM in the universe) and ii) dwarf spheroidal galaxies (the cleanest, nearby 
and bound DM halos). 



Dark Matter Probes 

Direct detection of DM particles is the cleanest and most decisive discriminant (see, 
e.g., Munoz 2003 for a review). However, it would be especially interesting if astronom- 
ical techniques were to reveal some of the fundamental particle properties predicted by 
fundamental theories. 

The dark side of the universe sends us, in fact, signals of the presence and of the na- 
ture of DM that can be recorded using different astrophysical probes. These probes are 
of inference and physical character. 

Inference probes [i.e., the CMB anisotropy spectrum (see, e.g., Hu & Dodelson 2002, 
Spergel et al. 2003, Komatsu et al. 2010), the dynamics of galaxies (Zwicky 1933), the 
hydrodynamics of the hot intra-cluster gas (see Sarazin 1988, Arnaud 2005 for a re- 
view) and the gravitational lensing distortion of background galaxies by the intervening 
potential wells of galaxy clusters (see Bartelmann & Schneider 1999 for a review and 
references therein)] tell us about the presence, the total amount and the spatial distri- 
bution of DM in the large scale structures of the universe but cannot provide detailed 
information on the nature of DM. 

Physical probes tell us about the nature and the physical properties of the DM particles 
and can be obtained by studying the astrophysical signals of their annihilation/decay in 
the atmospheres of DM-dominated structures (like galaxy cluster and galaxies). These 
probes can be recorded over a wide range of frequencies from radio to gamma-rays and 
prelude to a full multi-frequency, multi-experiment and multi-messenger search for the 
nature of DM in cosmic structures. 



A test case: neutralino DM 



Among the viable competitors for having a cosmologically relevant DM species, the 
leading candidate is the lightest particle of the minimal supersymmetric extension of the 
Standard Model (MSSM, see Jungman et al. 1996), plausibly the neutralino x, with a 
mass M x in the range between a few GeV to a several hundreds of GeV (see Baltz 2004 
for a review). Information on the nature and physical properties of the neutralino DM 
can be obtained by studying the astrophysical signals of their interaction/annihilation 
in the halos of cosmic structures (galaxy clusters and/or galaxies). These signals 5 
involve, in the case of a % DM, emission of gamma-rays, neutrinos, together with the 
synchrotron and bremsstrahlung radiation and the Inverse Compton Scattering (ICS) of 
the CMB (and other background) photons by the secondary electrons produced in the 
DM annihilation process (see Fig.l). 



FIGURE 1. A simple model which shows the basic astrophysical mechanisms underlying the search for 
the nature of % DM particles through the emission features occurring in large-scale structures (e.g., galaxy 
clusters and galaxies). These mechanisms are, among others: y-ray emission from n° — > y+ y, relativistic 
bremsstrahlung of secondary e ± and Inverse Compton Scattering (ICS) of CMB photons by secondary e ± ; 
X-ray /UV emission due to non-thermal bremsstrahlung and ICS of background photons by secondary e ± ; 
synchrotron emission by secondary e diffusing in the intra-cluster magnetic field; Sunyaev-Zel'dovich 
(SZom) effect (i.e. ICS of CMB photons by secondary e*) 

The x annihilation rate is R = n x (r)(ov), where n x (r) = n X fig(r) is the neutralino 
number density with radial distribution given by the function g(r), and (av) is the XX 
annihilation cross section averaged over a thermal velocity distribution at freeze-out 



Neutralino DM annihilation produces several types of particle and anti-particle fluxes, whose complete 
description is not discussed here for the sake of brevity. We refer the interested reader to Colafrancesco, 
Profumo & Ullio (2006) for the case of galaxy clusters and Colafrancesco, Profumo & Ullio (2007) for 
the case of dwarf galaxies. 



v 




temperature (Jungman et al. 1996). The range of neutralino masses and pair annihilation 
cross sections in the most general supersymmetric DM setup is extremely wide. Neu- 
tralinos as light as few GeV (see Bottino et al. 2003), and as heavy as hundreds of TeV 
(see Profumo 2005) can account for the observed CDM density through thermal pro- 
duction mechanisms, and essentially no constraints apply in the case of non-thermally 
produced neutralinos. Turning to the viable range of neutralino pair annihilation cross 
sections, coannihilation processes do not allow to set any lower bound, while on purely 
theoretical grounds a general upper limit on (<7v) ^ 10~ (M^./TeV) _2 cm 3 /s can be 
set (Profumo 2005). The only general argument which ties the relic abundance of a neu- 
tralino (WIMP) with its pair annihilation cross section is given by 

a DM h 2 ~ (3 x 10" 27 cm 3 /s)/((Jv) (1) 

(see Jungman et al. 1996). This relation can be, however, badly violated in the general 
MSSM, or even within minimal setups, such as the minimal supergravity scenario (see 
discussion in Colafrancesco, Profumo & Ullio 2006). 



Neutralino annihilations in DM halos 

Neutralinos which annihilate inside a DM halo produce quarks, leptons, vector bosons 
and Higgs bosons, depending on their mass and physical composition. Electrons and 
positrons (hereafter refereed to as electrons for simplicity) are then produced from the 
decay of the final heavy fermions and bosons. The different composition of the XX 
annihilation final state will in general affect the form of the electron spectrum. 

Neutral pions produced in %% annihilation decay rapidly in n° — > yy and generate 
most of the continuum spectrum at energies E ^ 1 GeV. 

Secondary electrons are produced through various prompt generation mechanisms 
and by the decay of charged pions, ->■ /i ± v M (v M ), with ju* ->• + v M (v M ) + v e (v e ) 
and produce e ± , muons and neutrinos. 

Secondary electrons are subject to spatial diffusion and energy losses. Both spatial diffu- 
sion and energy losses contribute to determine the evolution of the source spectrum into 
the equilibrium spectrum of these particles, i.e. the quantity which is used to determine 
the multi-frequency spectral energy distribution (SED) induced by DM annihilation. The 
time evolution of the secondary electron spectrum is described by the transport equation: 

= V [DVn e ] + — [b e (E)n e ] + Q e (E, r) , (2) 

where Q e (E,r) is the e^ source spectrum, n e (E,r) is the e ± equilibrium spec- 
trum and b e (given here in units of GeV/s) is the e ± energy loss per unit 
time b e = b IC s + b synch + b brem + b Cou l , with b 1C s ~ 2.5 • 10~ 17 '(E /GeV) 2 , 
b synch » 2.54- lO- l *B 2 (E/GeV) 2 , b brem w 1.51 • \Q- l \n th / cm^) (log(r/n, A ) +0.36), 

bcoul ~ 7 • 10~ 16 \n t hj 'cm' 3 ) (1 +log(r/n f ^)/75). Here n th is the ambient gas density 
and r = E/m e c 2 . 

The DM source spectrum, Q e (E, r), is constant over time, and under the assumption that 



the population of high-energy e ± can be described by a quasi- stationary {dn e /dt ~ 0) 
transport equation, the secondary electron spectrum n e (E,r) reaches its equilibrium 
configuration mainly due to synchrotron and ICS losses at energies E ^ 150 MeV and 
to Coulomb losses at smaller energies. 

The diffusion coefficient D in eq.(2) sets the amount of spatial diffusion for the sec- 
ondary electrons: it turns out that diffusion can be neglected in galaxy clusters while it is 
relevant on galactic and sub-galactic scales (see discussion in Colafrancesco, Profumo 
& Ullio 2006, 2007). 

To get a more physical insight on the relevance of spatial diffusion in large-scale struc- 
tures, it is useful to consider the following qualitative solution for the average electron 
density 

dnJE.r) _ . , V s Tn 

L « [Qe(E,r)T l0SS ] x — J— x O) 
dE V s + V Id + Tioss 

(see Colafrancesco 2005, Colafrancesco et al. 2006) which resumes the relevant features 
of the transport equation (eq. 2). Here, V s « R\ and V X 3 (E) are the volumes occu- 
pied by the DM source and the one occupied by the diffusing electrons which travel a 
distance X(E) m [D(E) ■ Xi 0SS (E)] 1 / 2 before loosing much of their initial energy. The rel- 
evant time scales in eq. (2) are the diffusion time-scale, Xd ~ RjjD(E), and the energy 
loss time-scale Xi oss = E/b e (E), where we assume the generic scaling of the diffusion 
coefficient D(E) = D (E/Eo) r -B" r - 

For E > E* = (DoEo/R^boB^) 1 ^ 1 ^ [for simplicity we kept leading terms only, im- 
plementing b(E) ~ bo(Bf 1 )(E /GeV) 2 + bcoul]^ the condition Xd > X[ oss (and consis- 
tently X(E) < Rfr) holds, the diffusion is not relevant and the solution of eq. (2) is 
dn e /dE ~ Q e {E : r)xi 0SS and shows an energy spectrum ~ Q(E) E ' . This situation 
(X(E) < Rh, Xd > Xioss) applies to the regime of galaxy clusters, i.e. structures on ~ 
Mpc scales (see Colafrancesco et al. 2006). 

For E < E*, the condition Xd < X[ oss (and consistently X(E) > Rh) holds, the diffusion 
is relevant and the solution of Eq. (2) is dn e /dE ~ [Q e {E,r)xD] x (V s /V ) and shows 
an energy spectrum ~ Q(E) ■ £'( 2_5 ^)/ 2 which is flatter or equal to the previous case for 
reasonable values 7= 1/3 — 1. This last situation (X(E) > Rh, Xd < Xi oss ) applies to 
the regime of dwarf galaxies, i.e. structures on ~ kpc scales (see Colafrancesco 2005, 
Colafrancesco et al. 2007). 

Secondary electrons eventually produce radiation by synchrotron in the magnetized 
atmosphere of cosmic structures, bremsstrahlung with ambient protons and ions, and 
Inverse Compton Scattering (ICS) of CMB (and other background) photons (and hence 
an SZ effect, Colafrancesco 2004). These secondary particles also produce heating of 
the ambient gas by Coulomb collisions with the ambient plasma particles. 



Spectral Energy Distribution from DM annihilation 

The astrophysical signals of neutralino DM annihilation computed in various DM 
models can be visible over the entire e.m. spectrum, from radio to y-ray frequencies 
(see Figs. 2 and 3). As pointed by Colafrancesco et al. (2006), the relevant physical 
properties which determine the features of the emitted radiation are the composition 



of the neutralino, its mass, and the value of the annihilation cross section. We consider 
here, for the sake of illustration, a few representative DM models with low (40 GeV), 
intermediate (81 GeV) and high (100 GeV) neutralino mass. 

Gamma rays. Gamma-ray emission is predominantly due to the hadronization of the 
decay products of XX annihilation with a continuum y-ray spectrum due to the decay 
7t° — > 7+ 7 (see, e.g., Colafrancesco & Mele 2001, Colafrancesco et al. 2006), even 
though the direct neutralino annihilation results in a line emission at an energy ~ M x . 
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FIGURE 2. Multi-frequency SED of the Coma cluster for two representative DM models: M x — 40 GeV 
(bb; left) andM x =81 GeV (W + W~; right). The halo profile is a NFW profile withM,, lr = 0.9 IO^Mq/T 1 
and c V ir = 10, with DM subhalo setup as given in Colafrancesco et al. (2006). The scaling of the multi- 
frequency SED with the value for the mean magnetic field in Coma is shown for the two DM models. 
The neutralino pair annihilation rate has been tuned to fit the radio halo data Figures from Colafrancesco 
et al. (2006). 

Gamma-ray emission is also expected from secondary e ± through bremsstrahlung and 
ICS of CMB photons up to high energies (see Figs. 2, 3). 

For the case of the Coma cluster, the gamma-ray flux produced by the low mass model 
(M x = 40 GeV, bb) is dominated by the continuum n° — > yy component and it is a factor 
~ 5 lower than the EGRET upper limit of Coma at its peak frequency once the annihi- 
lation cross section is normalized to fit the radio halo data (see Fig. 2, left panel). Such 
gamma-ray flux could be, nonetheless, detectable by the Fermi (GLAST)-LAT experi- 
ment for low values of the magnetic field ^0.5 /lG. A DM model with intermediate 
mass (M x = 40 GeV, W^) predicts lower gamma-ray flux below the EGRET limit but 
still detectable by Fermi (GLAST)-LAT only for very low values of the magnetic field 
^ 0.2 jiG (Fig 2, right panel). 

The rather low neutralino masses of these models make them difficult to be testable by 
Cherenkov gamma-ray detectors operating at high threshold energies. 

For the case of smaller cosmic structures, like the Draco dwarf galaxy, a DM model 
with M x = 100 GeV (normalized to the EGRET upper limit for Draco) predicts that 
the dominant gamma-ray emission is still given by the continuum n° — > yy component 
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FIGURE 3. Multi-frequency SED of Draco for a 100 GeV neutralino annihilating into bb and for a 
magnetic field of 1/iG (left). The effect of varying the magnetic field strength is also shown (right). The 
neutralino pair annihilation rate has been tuned to give a gamma-ray signal at the level of the EGRET 
upper limit. Figures from Colafrancesco et al. (2007). 



while the dominant ICS emission (i.e., ICS on CMB photons) is a factor ~ 10 2 less 
intense, a fact mainly due to diffusion effects. The other ICS component considered for 
Draco (i.e., ICS on starlight photons) peaks at frequencies comparable with those of the 
K° —7- 77 component, but is negligible being a factor ~ 10 4 less intense (see Fig. 3). 
Also this model with M x = 100 GeV is difficult to be tested by Cherenkov experiments. 

Radio emission. Secondary e ± produced by XX annihilation can generate synchrotron 
emission in the magnetized atmosphere of galaxy clusters (as well as galaxies) 
which could be observed at radio frequencies as a diffuse radio emission (i.e. a 
radio halo or haze) centered on the DM halo. Observations of cluster radio-halos 
are, in principle, very effective in constraining the neutralino mass and composition 
(Colafrancesco & Mele 2001, Colafrancesco et al. 2006), under the hypothesis that 
DM annihilation provides a major contribution to the radio-halo flux. Under this hy- 
pothesis, a pure energy requirement requires that the neutralino mass is bound to be 
M x > 23.4GeV(v/GHz) 1 / 2 (5/^G) -1 / 2 in order that the secondary e ± emit at frequen- 
cies V > 1 GHz, as observed in cluster radio halos (see Fig. 2). A soft DM model (bb 
with M x = 40 GeV) is able to reproduce both the overall radio-halo spectrum of Coma 
and the spatial distribution of its surface brightness (see Fig.4), while a hard DM model 
with M x = 81 GeV) is excluded being its radio spectrum too flat to reproduce 
the Coma data. 

For the case of dwarf galaxies (e.g. Draco), radio emission is strongly affected by 
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FIGURE 4. Left. The radio flux density spectrum of Coma in two DM models: soft spectrum due 
to a bb annihilation final state (solid line) with M x = 40 GeV) and hard spectrum due to a 
channel (dashed line) with M x —81 GeV. Right. Surface brightness distribution of Coma at v = 1 .4 GHz, 
within a beam equal to 9'.35 (HPBW), for the the soft neutralino model (solid). We consider a magnetic 
field B(r) = B {1 + (r/r cl ) 2 ) • [1 + (r/r t2 ) 2 ]^ with B = 0.55 pG, J3 = 2.7, r cl = 3' r c2 = 17.'5. The 
contributions to the radio brightness from the smooth DM halo component only (dotted) and from subhalos 
only (dashed) are also displayed (see Colafrancesco et al. 2006 for details). 



propagation effects. 6 Fig. 5 shows that for a propagation set up (set up #1) with a 

Kolmogorov spectrum (D E 1 / 3 ^^ 1 ) there is a depletion of the electron populations 
with a significant fraction leaving the diffusion region, while for a propagation set up 
(set up #2) with a steeper spectrum (D E ( 'BJ L ( ') they are more efficiently confined 
within the diffusion region but still significantly misplaced with respect to the emission 
region. As a consequence, also the spectral shape of the radio flux of Draco is affected by 
diffusion effects which produce a steeper spectral slope when the electron populations 
are more efficiently confined within the diffusion region (i.e. propagation set up #2) with 
respect to the case (i.e. propagation set up #1) where there is a depletion of the electron 
populations with a significant fraction leaving the diffusion region (see Colafrancesco et 
al. 2007 for details). 



6 For the diffusion coefficient we consider here the case of a Kolmogorov form D(E) = 

Do/B 1 / 3 (E/l GeV) 1/3 , with D Q = 3 • 10 28 cm 2 s _1 , in analogy with its value for the Milky Way (here 
Bp is the magnetic field in units of flG). The dimension of the diffusion zone is, consistently with the 
Milky Way picture, about twice the radial size of the luminous component, i.e. ps 102 arcmin for Draco 
(set #1 of propagation parameters). An extreme diffusion model in which the diffusion coefficient is de- 
creased by two orders of magnitudes down to Do = 3 • 10 26 cm 2 s _1 (implying a much smaller scale of 
uniformity for the magnetic field), and with a steeper scaling in energy, D(E) — Dq (E/l GeV) -0 6 (this is 
the form sometime assumed for the Milky Way) is considered for comparison (we label this propagation 
parameter configuration set #2) 
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FIGURE 5. Radio flux density spectrum (left) and surface brightness distribution (right) of Draco at 
V = 4.9 GHz for a sample % model with M x = 100 GeV annihilating into bb with a cross section tuned 
as to give 2 events in the EGRET experiment. Results are shown for the two choices of propagation 
parameters, as labelled; the case when spatial diffusion is neglected is also shown for comparison in the 
right panel (dotted curve). Surface brightness radial profiles are plotted for a 3 arcsec angular acceptance, 
corresponding to the VLA resolution used to search for a point-like radio source at the center of Draco 
(the VLA upper limit is also plotted in the figure). 



ICS of CMB: from microwaves to gamma-rays Secondary e produced by neutralino 
annihilation up- scatter CMB (and other background) photons that will redistribute over 
a wide frequency range up to gamma-ray frequencies. A soft neutralino model with 
M x = 40 GeV and (ov) = 4.7 ■ 1(T 25 cm 3 * -1 , with = 1.2 yields EUV and HXR 
fluxes which are more than one order of magnitude fainter than the Coma data (see 
Fig. 2). Increasing M x does not provide a decent fit of the radio-halo spectrum (see 
Fig. 2, right panel) and yields, in addition, extremely faint EUV, HXR and gamma- 
ray fluxes, which turn out to be undetectable even by Fermi (GLAST) and/or by the 
next coming high-energy experiments. Increasing (<7v) by a factor ~ 10 2 (i.e., up to val- 
ues (ov) ~ 7 • 10 _23 cm 3 5 _1 ) can fit the EUV and HXR data on Coma but the relative 
71° — > yy gamma-ray flux exceeds the EGRET upper limit on Coma. 
Lowering the magnetic fields down to values ~ 0.15/iG can fit both the HXR and the 
EUV fluxes of Coma but also in this case the n° — > YY gamma-ray flux predicted by the 
same model exceeds the EGRET limit on Coma, rendering untenable also this alterna- 
tive. Actually, the EGRET upper limit on Coma set a strong constraint on the combina- 
tion of values B and (ov) so that magnetic field larger than ^ 0.3 /lG are required for 
the parameter setup of the bb model with M x = 40 GeV (Fig. 22 of Colafrancesco et al. 
(2006) shows the upper limits on the value of (ov) as a function of the assumed value 
of the mean magnetic field of Coma). 

According to these results, it is impossible to fit all the available data on Coma for a 
consistent choice of the DM model and of the cluster magnetic field (see also discussion 
in Colafrancesco et al. 2010). 



For the Draco dwarf galaxy the dominant ICS on CMB component produces fluxes 
of X-rays at the level of ~ 10~ 15 — 10~ 14 erg cm~ 2 s _1 when the gamma-ray flux is 
normalized to the EGRET upper limit. The constraints obtainable by Fermi (GLAST) 
observation of Draco will hence set much more realistic expectations for the diffuse X- 
ray emission produced from DM annihilation in Draco which could eventually be tested 
with high sensitivity X-ray, and especially HXR and soft y-ray observations where this 
ICS spectrum peaks at E ~ 100 keV. 

ICS of CMB: SZ effect from DM annihilation. Secondary e ± produced by neutralino 
annihilation interact with the CMB photons and up- scatter them to higher frequencies 
producing a peculiar SZ effect (as originally realized by Colafrancesco 2004) with 
specific spectral and spatial features. The DM induced spectral distortion writes as 
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A/dm to = 2^J-vdm g(x) , (4) 



where Tq is the CMB temperature and the Comptonization parameter is 

Ot f 

vdm = j / Poudi , (5) 

m e c z J 

in terms of the pressure i^M contributed by the secondary e ± . The function g(x), with 
x = hv/k^TQ, for the DM produced secondary e ± can be written as 



s( x ) = -T^{t— I io( xe s )P(s)ds-i (x) 



(6) 



in terms of the secondary e ± optical depth Xdm = &t J din e (with n e given by the 
solution of eq.2), of the photon redistribution function Pis) = J dpf e (p)P s (s; p) with s = 
ln(v'/v), in terms of the photon frequency increase factor v'/v, of the e ± momentum 
(p) distribution f e (p), of io(x) = 2(fcB?b) 3 / (he) 2 -x 3 / (e x — 1), and of the quantity 

{e) = ^-f Pdi = = r dpU{p)\pv{p)m t c (7) 

Tdm J J n s dl Jo 3 

which is the average energy of the secondary electron population (see Colafrancesco 
2004). Fig. 6 shows the maps of the CMB temperature change, 
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at 223 GHz as produced by the SZom effect in the neutralino models with mass of 40 
GeV (bb) and 81 GeV (VF ± ), compared to the temperature change due to the residual 
thermal SZ effect produced by the intracluster gas of the bullet cluster. The spatial 
separation between the DM and thermal SZE of this cluster is due to the fact that 
the SZpM effect has a very different spectral shape with respect to the thermal SZ 
effect: it yields a temperature decrement at all the microwave frequencies ^ 600 GHz, 
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FIGURE 6. The simulated SZ maps of the cluster 1ES0657-556 observable at v = 223 GHz with a 
telescope of angular resolution similar to SPT (see Colafrancesco et al. 2007 for details) for a neutralino 
mass M x = 40 GeV, bb (left panel) and M x = 81 GeV W ± (right panel) with the DM density model setup 
of Colafrancesco et al. (2006). 



where the thermal SZ effect is predominantly observed, and produces a temperature 
increase only at very high frequencies > 600 GHz. This behavior is produced by the 
large frequency shift of CMB photons induced by the ICS off the relativistic secondary 
electrons generated by the neutralino annihilation. As a consequence, the zero of the 
SZom effect is effectively removed from the microwave range and shifted to a quite high 
frequency ~ 600 GHz with respect to the zero of the thermal SZ effect, a result which 
allows, in principle, to estimate directly the pressure of the two electron populations and 
hence to derive constraints on the neutralino DM model (see Colafrancesco 2004). 
It is, however, necessary to stress that in such frequency range there are other possible 
contributions to the overall SZ effect, like the kinematic SZ effect and other non-thermal 
SZ effects which could provide additional biases. 

We stress that a SZdm effect is expected in every neutralino DM halo and its amplitude 
depends basically on the optical depth of the DM-produced secondary electrons, and 
hence on the detailed distribution of the equilibrium spectrum of the secondary electrons. 
The SZdm effect produced in a dwarf galaxy, like Draco (see, e.g., Colafrancesco 2004, 
Culverhouse, Ewans and Colafrancesco 2006), is however quite low when the spatial 
diffusion of secondary electrons is efficiently operating: we find, in fact, that the SZ 
signal towards the center of Draco is negligible even when we normalize the gamma-ray 
signal at the level of the EGRET upper limit. 

Heating. Low energy secondary electrons produced by neutralino annihilation heat the 
intra-cluster gas by Coulomb collisions since the Coulomb loss term dominates the 
energy losses at E ^ 150 MeV. The DM-induced heating rate at the center of galaxy 
clusters is usually higher than the intra-cluster gas cooling rate for steep DM density 
profiles (like the NFW one). The radius of the region in which such steep DM density 
profile produce an excess heating increases with increasing neutralino mass, once the 
DM-produced SED are normalized to fit the radio halo data. The heating effect provides, 
hence, strong constraint to the annihilation cross-section for neutralino DM models. In 



order to have the DM-induced heating rate lower than the bremsstrahlung cooling rate 
at the cluster center the annihilation cross section must be reduced by a large factor 
(see Colafrancesco et al. 2010 for a discussion in the case of Coma cluster). Cored 
DM density profiles alleviate the over-heating problem for neutralino DM models by 
reducing the electron equilibrium spectrum n e (E,r) in the cluster central region. This 
makes DM models more consistent with the cluster heating rate constraints for low 
(a few to ~ 10 GeV) and intermediate (~ 40 — 60 GeV) neutralino mass, but not 
for high neutralino mass (~ 500 GeV or more) that still produce excess heating (see 
Colafrancesco et al. 2010). 

Multi-Messengers. Stronger constraints to DM annihilation models can be set by a 
multi-messenger analysis of DM signals. Such multi-messenger constraints include, e.g., 
observations and limits on positrons, antiprotons, radio and y-rays from the Galactic 
Centre region and the optical depth of CMB photons (see discussion in Colafrancesco 
et al. 2010); multi-frequency (radio, microwave, X-ray and gamma-ray) observations 
from satellites of the MW, external galaxies (like M31, e.g. Bordello et al. 2009) and 
cosmic background radiation, in addition to those obtained from dwarf galaxies and 
galaxy clusters (see Colafrancesco 2006, Profumo & Ullio 2010 for reviews). 
We stress that these multi-messenger bounds are quite robust and not easily avoidable, 
and therefore they set further constraints to any specific multi-frequency analysis of DM 
models, as we have described in this review. 

Cosmic rays from neutralino annihilation 

Neutralino annihilation in nearby DM clumps (like e.g. the Galactic center or galactic 
DM satellites) produces cosmic rays that diffuse away and can be directly recorded by 
cosmic rays experiments. 

Some of the cosmic ray observations (at energy between 10 GeV and a few TeV) ob- 
tained with PAMELA (Adriani et al. 2009), ATIC (Chang et al. 2008), Fermi (Abdo et 
al. 2009), and HESS (Aharonian et al. 2009) show a positron excesses over background 
expectations (Strong et al. 2009). 

An attractive and widely explored possibility is that the observed positron excesses 
are produced by neutralino (WIMP) DM annihilation, even though these excesses have 
rather plausible astrophysical explanations (see, e.g., Hooper et al. 2009, Yuksel et al. 
2009, Profumo 2008, Dado & Dar 2009, Biermann et al. 2009, Katz et al. 2009). How- 
ever, if the neutralino annihilation cross section is of the order of thermally-averaged 
annihilation cross section at freeze out, i.e. (ctv) ^3x 10~ 26 cm 3 /s, the resulting anni- 
hilation signal is far too small to explain the observed cosmic ray excesses. 
A seemingly attractive solution is to postulate that DM interacts with a light force carrier 
with fine structure constant ax = A 2 / (4k) (see, e.g., Cirelli et al. 2009). This effec- 
tively enhances the annihilation cross section by a quite large factor S = i_ e ™-xax/v t ) 
(usually referred to as Sommerfeld enhancement, see e.g. Arkani-Hamed et al. 2009), 
where v re i is the DM particle relative velocity. Since the required DM mass to explain 
the positron observations is ~ TeV, the required is ~ GeV, values S ~ 10 3 can be 
obtained, assuming (ov re i) m (ctv), that can fit the positron excess. The Sommerfeld 
enhancement therefore provides an elegant mechanism for boosting neutralino DM an- 



nihilations now. 

Of course, for a viable solution, DM must not only annihilate with the correct rate, but 
it must also be produced with the right density and form all large-scale structure in 
agreement with observations. However, Feng et al. (2009) have shown that the required 
enhancement implies thermal relic densities that are too small to be all of DM. In ad- 
dition, upper bounds on possible Sommerfeld enhancements can be derived from the 
observation of elliptical galactic dark matter halos (like the case of NGC270), and these 
bounds also (generically) seem to exclude enhancements that can explain the observed 
positron excesses (see discussion by Feng et al. 2009). 



Other DM scenarios 

Recent progress has greatly expanded the list of well-motivated candidates and the 
possible signatures of DM (see Feng 2010 for a recent review). Beyond any attempt to 
be exhaustive and/or complete in this respect, we briefly discuss here the interesting case 
of a light DM candidate, i.e. sterile neutrinos. 

Sterile neutrinos may be produced in a number of ways and their relic density 
depends on the sterile neutrino mass and mixing angle, but all of the mechanisms require 
small masses m s and mixing angles for sterile neutrinos to be a viable DM candidate. 
As for the astrophysical (indirect detection) search, it is interesting to notice that the 
radiative decay of sterile neutrinos, v s — >■ V; + y (where v, indicate the standard low-mass 
neutrinos), produces a narrow line emission whose energy provides information on the 
sterile neutrino mass m s . Therefore, X-ray spectral observations from galaxy clusters 
are a powerful tool to set contraints on sterile neutrinos in the plane m s — sin 2 (26). 
The available constraints on sterile neutrinos from X-ray spectra of galaxy clusters, 
combined with those obtained from the Cosmic X-ray Background, Lya limits and 
gamma-ray line limits from the MW are shown in Fig. 7. The constraints from Coma 
observations in the 20-80 keV band are shown by the cyan dashed area. The possible 
interpretation of the intensity excess in the 8.7 keV line (at the energy of the FeXXVI 
Ly-y line) in the spectrum of the Galactic center observed by the Suzaku X-ray mission 
in terms of decay of a sterile neutrino with mass of 17.4 keV and value of the mixing 
angle sin 2 (29) = (4.1 ±2.2) • 10~ 12 (see Prokhorov & Silk 2010), lies in the allowed 
region of mass-mixing angle for DM sterile neutrino shown is Fig. 7. 
This figure shows that models with lower mixing angles and neutrino masses m s up 
to a few hundreds keV or even ^ MeV are still available. In this case, next generation 
high-sensitivity hard X-ray detectors like NuSTAR and/or NeXT, or next coming soft 
gamma-ray experiments will be able to set relevant constraints to this light DM model. 

DARK MATTER IN MODERN COSMOLOGY: THE FUTURE 

All solid evidence for DM is gravitational, and there is also strong evidence against DM 
having strong or electromagnetic interactions. However, direct and indirect probes for 
DM have, so far, not yet given a definite (positive) answer. In addition, some of the 
puzzling anomalies (e.g., DAMA, PAMELA) are not easy to explain within canonical 
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FIGURE 7. The sterile neutrino mass m s and mixing sin 2 (29) parameter space, where the shaded 
regions are excluded. The strongest direct bounds are shown, labeled as Milky Way (Yuksel et al. 2008), 
Cosmic X-ray Background (Boyarsky et al. 2006), and X-ray limits (summarized by Watson et al. 2006). 
The strongest indirect bounds (Seljak et al. 2005), Viel et al. 2006) are shown by the grey horizontal band. 
The excluded Dodelson-Widrow model (Dodelson & Widrow 1994) is shown by the solid line; rightward, 
the DM density is too high (stripes). The dotted lines indicate example models (see Abazajian & Fuller 
2002), now truncated by the available constraints. The cyan shaded area indicates our constraints from the 
HXR (soft gamma-ray) limit on Coma (Colafrancesco 2008). Figure from Colafrancesco (2008) adapted 
from the original one by Yuksel et al. (2008). 



WIMP DM models (see Feng 2009 for a discussion). This experimental frustration and 
theoretical embarassement have motivated the search for new DM candidates. Among 
the logical alternatives to the WIMP paradigm (i.e., the fact that particle physics theories 
designed to explain the origin of the weak scale often naturally contain particles with the 
right relic density to be DM), then, one of the most widely explored ones is provided by 
hidden DM, that is, DM that has no standard model gauge interactions (see Feng 2009, 
2010 for reviews). 

At the moment, therefore, the DM search is continuing to be a great challenge. The 
direction in this effort are: direct search experiments, indirect probes of DM signals 
coming from cosmic sources, indirect indications on the nature of DM coming from 
laboratory experiments, like the LHC experiment. The next future is certainly bright in 
this context and expectations are rather high. 

In the cosmological and astrophysical context, the refinement for the indirect DM 
search calls for a Multi 3 approach in the optimal astrophysical laboratories: i) multi- 
frequency to probe the consistency of DM annihilation/decay signals across the e.m. 
spectrum; ii) multi-messenger to cross-check the possible e.m. signals in cosmic struc- 
tures at different scales; iii) multi-experiment to determine the robustness of the DM 
signals with different experimental techniques. 



Multi 3 DM search: optimal astrophysical laboratories 



The analysis of the spatial and spectral intensity of the astrophysical signals coming 
from neutralino annihilation is a powerful tool to unveil the elusive nature of DM. 
However, the DM-induced signals are expected to be confused or even overcome by 
other astrophysical signals originating from the ambient gas and/or from the relativistic 
plasmas present in the atmospheres of galaxy clusters and galaxies, especially when all 
these components are co-spatially distributed with the DM component. 
An ideal system to detect DM annihilation signals would be a system which is either 
devoid of diffuse emitting material (this is the case of dark galaxies, like many dwarf 
spheroidals) or a system with a clear spatial separation between the various matter 
components (this is, indeed, the case of the cluster 1ES0657-556 where the spatial 
distribution of DM which is clearly offset with respect to that of the intracluster gas). 

A multi-frequency analysis greatly helps to disentangle DM signals in these optimal 
laboratories from signals of different astrophysical origin. 

In the multi-messenger analysis previously evoked, dark (dwarf) galaxies are among 
the best sites for the astrophysical search for the nature of DM but the relative multi- 
frequency SED are usually quite dim. Definite probes of DM signals in such systems 
have to be then complemented with probes coming from other DM halos like galaxy 
clusters (on larger scales) and the center of the Galaxy (on closer scales). 

The multi-experiment strategy combining radio and gamma-ray observations of 
dwarf galaxies, like Draco, obtainable with high sensitivity instruments (SKA, LOFAR, 
EVLA, Fermi) could set strong constraints on the nature of the DM particles see 
Colafrancesco et al. 2007 for a discussion). Fig. 8 shows the limits in the M x — (<7v) 
plane set from a Multi 3 analysis of neutralino annihilation signals (taking into account 
the Draco dwarf galaxy and some of the signals coming from the Milky Way) where the 
most constraining frequency ranges are, in fact, the radio and the gamma-ray bands. 

As for galaxy clusters, a Multi 3 approach indicates that in the case of 1ES0657-556 the 
expected gamma-ray emission associated to the DM clumps is too low ( ^ 1 count vs. 
~ 10 background counts at E > 1 GeV) and cannot be resolved by Fermi from other 
possible sources of gamma-ray emission, both from the cluster 1ES0657-556 and from 
AGNs (or other y-ray emitting galaxies) in the field. Radio telescopes have, in principle, 
excellent resolution and sensitivity to probe the different spectra and brightness distri- 
bution of the DM-induced synchrotron emission, but the theoretical uncertainties asso- 
ciated to the radio emission from the DM clumps in 1ES0657-556 render the interpre- 
tation of the expected signals quite uncertain: we evaluated (Colafrancesco et al. 2007) 
that the DM induced synchrotron emission from the largest DM clump is ~ 3 — 10 mJy 
(for a smooth or smooth plus 50% mass dumpiness NFW DM profile, soft DM model 
with M x = 40 GeV with a B = 1 /iG) at v = 100 MHz, still marginally detectable by 
LOFAR. In such a context, the possible detection of the SZom effect from this system, 
with the next generation high- sensitivity and high-resolution experiments, will provide 
an important complementary, and maybe unique, probe of the nature of DM. 
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FIGURE 8. The limits in the M x — (pv) plane set by current and future gamma-ray experiments in 
a Multi 3 analysis of a bb neutralino annihilation model. Left. The solid black line indicates values of 
(<7v) required for a gamma-ray flux of the Draco dwarf galaxy matching the 2 events in EGRET between 
1 and 10 GeV, for exposures and angular cuts as specified in the data analysis, assuming a reference 
NFW halo model. The black dotted and dot-dashed curves show the limits from the flux of, respectively, 
antiprotons and positrons in the Milky Way halo (the two lines corresponding to two different cosmic ray 
propagation setups); The solid red, dark blue and light blue curves indicate the projected sensitivity of 
Fermi (GLAST), MAGIC and VERITAS for Draco. Right. The projected sensitivity of future diffuse 
radio source searches from the direction of Draco. For comparison, we also indicate the sensitivity 
of gamma-ray search experiments, and the constraints from the antiproton and positron flux (the two 
lines corresponding to two different propagation setups for the Milky Way). (Figures and details from 
Colafrancesco et al. 2007). 



DM or modified gravity? 

The DM particle solution of the global DM problem is not univocally accepted. 
A different, more radical approach to explain the cosmological DM problem can be taken 
if one notes that the evidence for missing mass arises because of a mismatch between the 
gravitational field one would predict from the observed mass distribution in the universe 
and the observed gravitational field. The observed discrepancies arise when the effective 
gravitational acceleration is around, or below, the value ciq ~ 10~ 7 cms~ 2 , that is in a 
regime of very weak gravitational field (see Ferreira & Starkman 2009 for a discussion). 
The point where modified gravity scenarios are at stake, is that the Newtonian theory of 
gravity - and general relativity - break down in this regime (see Fig. 9). 

It has been shown that modified theory of gravity that incorporates quantum effects 
can explain a crucial set of puzzling astronomical observations (galaxy flat rotation 
curves, galaxy cluster structure and the accelerated expansion of the universe), including 
the "wayward" motion of the Pioneer spacecraft in our solar system (e.g. Moffat 2005, 
2006; see also Brownstein & Moffat 2006a,b,c). 

Beyond the details of the specific approaches to a modified or extended theory of gravity 
that have been presented so far, it is important to notice that a relativistic theory of 



Solar System 



Globular Clusters 



Modified Gravity (DM) GR (no DM) 



Dwarf Galaxies 
Spiral Galaxies 



;Clusters of galaxies 



Large Scale Structure 

i i Cosmic Eixpansion 

1 I I II i i ' iTiiT | I i ii | 

10-w 10- 9 10- 9 10- 7 10- 3 10- 5 

Acceleration (in cm s~ a ) 

V 

FIGURE 9. Evidence for deviations from general relativity (GR) or evidence for the presence of Dark 
Matter (DM) tend to appear in systems in which the acceleration scale is weak (the left side, gray shaded 
part of the plot) at 7 x 10~ 8 cms~ 2 . There is strong evidence for either of the above in dwarf galaxies, 
spiral galaxies, galaxy clusters, large scale structures and in the expansion of the universe itself. Figure 
adapted from the original one by Ferreira & Starkman (2009). 



modified gravity allows to make a number of specific predictions on a wide range of 
scales, from the scales of compact objects such as stars or black holes (in the limit of 
strong gravitational field), to the scales of formed structures on large scales (such as 
galaxies and galaxy clusters) up to the very large-scale structure scales and the largest 
scale probed by CMB anisotropics (in the limit of weak gravitational field) (see, e.g., 
Ferreira & Starkman 2009, Capozziello et al. 2009, Capozziello & Francaviglia 2008, 
Schmidt et al. 2009, among others). 

It is, therefore, conceivable to use the next coming experimental probes on both strong 
gravitational fields (in the vicinity of very compact objects) and weak gravitational field 
(on large-scale cosmic structures) to set relevant constraints on the level and on the 
spatial scales at which the modification of the relativistic theory of gravity could occur. 

EPILOGUE 

Viable DM models which are consistent with WMAP limits on Q.dm and Q-de and with 
the structure and evolution of galaxy clusters are able to produce substantial astrophys- 
ical signals especially detectable at radio, microwave, X-ray and gamma-ray frequen- 
cies. The constraints that the multi-frequency observations of the optimal astrophysical 
laboratories for DM search can set on the (ov)-M x plane, combined with the WMAP 
(relic abundance) constraints on (<7v), are able to efficiently restrict the available neu- 
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tralino DM models. Additional restrictions of this plane may be obtained through a 
multi-messenger approach by comparing the previous astrophysical constraints to the 
constraints coming from both accelerator physics and from other cosmological probes 
(e.g., the study of the emission features from the galactic center region, galactic satel- 
lites, dwarf galaxies, external galaxies) which are sensitive to the amount and nature 
of the DM. Direct DM detection experiments have already explored large regions of 
the most optimistic SUSY models, and the planned next-generation experiments will 
probably be able to explore also the core of the SUSY models. In this context, the astro- 
physical study of DM annihilation proves to be complementary, and hardly competitive, 
especially when a full multi-messenger and multi-experiment approach is chosen. 
When combined with future accelerator results, such Multi 3 (multi-frequency, multi- 
messenger, multi-experiment) astrophysical search might greatly help us to unveil the 
elusive nature of Dark Matter. 

The extended theory of gravity, on both an empirical and theoretical level, seems to 
be the most viable alternative to a dark cosmological model in which Dark Matter will 
remain still elusive to laboratory and astrophysical probes. 
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